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Abstract 


Measurements have been made of the rising velocities and the shapes 
of spherical cap bubbles in a column 15cm in diameter. The liquids used 
were water, aqueous solutions of polyvinyl alcohol (PVA) with viscosities 
up to 216cP and ethyl alcohol; the gases were air, nitrogen and carbon 
dioxide. With carbon dioxide mass transfer coefficients were also obtained. 
Relationships between the measured quantities and the properties of the 
various liquids are considered. Liquid viscosities up to 50 cP have no effect on 
either the rising velocity or the shape of spherical cap bubbles. The addition 
of 0.5 % PVA to water results in a sharp drop in the mass transfer coefficient. 
Further additions cause it to increase and almost return to its original value. 


Messungen der Aufstiegsgeschwindigkeit und der Form von runden Blasen 
in einer Saule von 150mm innerem Durchmesser wurden durchgefiihrt. Die 
hierbei verwendeten Flissigkeiten waren Wasser, wasserige Losungen von 
Polyvinylalkohol (PVA) mit Viskosititen bis zu 216cP und Athylalkohol; 
die Gase waren Luft, Stickstoff und Kohlenséure. Bei Kohlenséure wurden 
auch Massentibertragungskoeffizienten ermittelt. Die Beziehungen zwischen 
den gemessenen Mengen und den Eigenschaften der diversen Flissigkeiten 
wurden berticksichtigt. Viskositaéten bis zu 50cP wirken sich weder auf die 
Aufstiegsgeschwindigkeit noch auf die Form der Kugelblasen aus. Ein Zusatz 
von 0.5% PVA zum Wasser ergibt scharfes Absinken des Masseniibertra- 
gungskoeffizienten. Weitere Zusatze haben neuerliches Ansteigen bis nahezu 
auf den urspriinglichen Wert zur Folge. 


2. MATERIALS 


1. INTRODUCTION 


The Importance of interactions between gas bubbles and liquids in many 
industrial processes is well known. Nevertheless, until recently, most of 
the investigations have been confined to small bubbles, less than lem’, 
Moreover, the liquid phases studied have generally been water or dilute 
aqueous solutions with similar properties. 

The investigation described in this paper was concentrated on the behaviour 
of large single bubbles of spherical cap shape rising in liquids of widely 
different properties. Since mass transfer between the bubble and the liquid 
may be influenced both by the rising velocity and the bubble shape, the initial 
part of the investigation was confined to these two factors. Subsequently, 
the mass transfer phenomena themselves were studied. 


2. MATERIALS 


Liquids employed in the investigation were tap water, distilled water, indus- 
trial ethyl alcohol and aqueous solutions of polyvinyl alcohol, made by the 
addition of ” Alcotex 88-40” powder to warm distilled water. The properties 
of these liquids are given in Table 1, where the references indicate literature 
values. The properties of the polyvinyl alcohol solutions were measured 
by standard techniques apart from the carbon dioxide diffusion coefficients 
which were determined by the laminar film absorption method developed by 
DAVIDSON and CULLEN, [8]. 

Bubbles were formed from cylinder air, nitrogen and carbon dioxide, the 
last two having less than 50 parts per million residual impurities, apart from 
moisture. 


Density 


Liquid eC ins?) 


Tap water (16 °C) 1.00 


Distilled water 1.00 
(16°C) 

Ethyl alcohol 0.79 
(17°C) (Industrial 
absolute) 


Aqueous polyvinyl 
alcohol (15°C) 
Percent PVA0.5% — 1.0! 


2.85 % 1.0! 
3.9% 1.0! 
5.4% 1.0! 
6.1% 1.0! 
' This work 


Surface Tension Viscosity 
(dyncm™*) (cP) 

73 1.11 

es 1.11 

22 1.25 

49.7 | on 

47.5 | 20° (12.5 °C) 
a7 0 43! 

46.9 | 115! 


a6 2-" 216" 47°C) 


Solubility of 


CO»x10° 


(mol cm °) 


4.39 
4.37+0.04! 


12.6[17| 


4.37+0.04" 
4.15" 
4.10" 
4.06" 
4.04" 


Diffusivity of 
CO2x10° 
(cm? s"*) 

1.37 (15°C) [8] 
1.3540.1 (15°C)! 


2.7 (16°C) [18] 


i3! ro) 
1.45' (15°C) 
1.4" (15°C) 
1.2' (15.5°C) 


0.15 


Table 1: Liquids employed during the investigation. (Except where stated, all properties from [13]) 


3. APPARATUS AND PROCEDURES 


3. APPARATUS AND PROCEDURES 


The bubbles were produced in 15cm i.d. Perspex tubes of lengths up to 
300 cm. Single bubbles were formed by the injection of gas into an inverted 
cup of spherical cap shape followed by manual rotation of the cup into the 
upright position. The larger bubbles showed a tendency to break up but, 
with care, this manual technique could be used to obtain perfect bubbles of 
volumes up to 50 cm”. 

Bubble volumes were determined by measuring the amounts of gas displaced 
from above the liquid column into a soap film meter, during the filling of 
the cup. Tests, in which the volume and pressure of the ingoing gas and the 
volume displaced above the column were measured, showed that the bubble 
volumes were obtained with a standard deviation of +0.02cm". The mean 
bubble velocities were determined by means of two photoelectric cells placed 
against the column, at different levels. The cells actuated a stopwatch via a 
solenoid so that transit times could be measured to within 0.01s. 

The shapes of bubbles rising in a 15cm square perspex tank were recorded 
photographically using an electronic flash, with a peak period of aa s. Mea- 
surements were taken from projections of the 35mm negatives. Comparisons 
between the volumes calculated from the linear dimensions so obtained and 
those measured directly indicated that the linear dimensions were accurate 
to approximately +4%. 

The rates of transfer of COg between bubbles and various liquids were deter- 
mined by measuring the changes in bubble volume which occurred during 
passage up the column. The experiments covered both absorption and des- 
orption. In absorption studies a predetermined volume of CO, was injected 
into the inverted cup and the displaced volume was measured precisely. The 
bubble was then released and the time of rise and the bubble volume at the 
surface were measured. Prior to each experiment both the liquid column 
and the gas space above were purged with nitrogen to ensure the system 
was free from COg,. In the desorption studies the column was saturated 
with CO, prior to each measurement by injecting a stream of small CO, 
bubbles (approx 2mm dia.) at the bottom. Nitrogen was then injected into 
the cup as quickly as possible in order to minimize the pickup of CO, prior 
to the bubble release. The time of rise and the initial and final volumes 
were measured as for the absorption studies. Fortunately a sudden change 
in the rate of movement of the soap film in the flowmeter occurred when 
the bubble broke surface. This facilitated precise measurement of the CO, 
transferred and the volume of the emergent bubble was taken to be that 
given by the flowmeter reading just prior to the change in rate. 

In the experiments using ethyl alcohol, all gases were saturated with alcohol 
vapour to reduce the effect of vapour transfer upon bubble volume. 


4. RISING VELOCITY 
4.1. Results 


Results of the velocity measurements are shown in Figure 1 where the equiva- 
lent radius is the radius of a sphere with the same volume as the spherical cap 
bubble. Since the change in hydrostatic pressure during the rise causes the 
bubble to expand, the particular volume chosen for calculating the equivalent 
radius was the arithmetic mean of the top and bottom values. 

The spread of results for the velocities of bubbles of identical size rising in 
water was found to be +1% of the mean value. Part of the deviation from 
the mean may be attributable to some bubbles not rising completely axially 
in the column. 

As predicted by the expressions of RIPPIN, [16] and of DAVIES and TayLor, 
[10], there was no measurable effect of liquid density. Though the experi- 
ments described here only covered the range from 0.79 gcm™® (ethyl alcohol) 
to 1 ecm” (water), in a further investigation DAVENPORT et al., [6], using 
mercury have shown that even a density of 13.5¢ cm’ does not affect the 
rising velocity. The mean terminal velocity is unaffected by surface tension 
in the range of 22dyncm™ (ethyl alcohol) to 73dyn cm | (water). The 
further investigation by DAVENPORT et al. referred to above showed that 
identical terminal velocities were obtained in mercury with a surface tension 
of 487 dynem™'. 

The results also show that surface active agents, namely small amounts 
(0.5 %) of polyvinyl alcohol, have no influence on the rising velocity of the 
bubbles. This confirms the results of previous workers, [2, 11] and it would 
appear such substances only influence smaller bubbles. 

Additions of PVA up to 4.5%, which increase the viscosities of the aqueous 
solutions to 54cP, have no appreciable effect on the rising velocity. Further 
additions to 115 cP cause some reduction in the velocity and this is even more 
marked when the viscosity is increased to 216cP. This effect is especially 
noticeable with bubbles of an equivalent radius less than 1 cm. Since the PVA 
behaves like a NEWTONian fluid, the effect of viscosity cannot be attributed 
to non-NEWTONian behaviour. The column length showed no appreciable 
influence on mean velocity. 


Velocity measurements on bubbles in water, which were either absorbing 
or desorbing CO, showed that the terminal velocity was not affected by mass 
transfer. 


4. RISING VELOCITY 


Curves displaced by 
minus 4 cm/s 


Rising velocity, cm/s 


ye) 12 
Radius of equivalent sphere, cm 


Figure 1: Bubble velocities vs. radius of equivalent sphere. 
Air in water in 70cm column +; air in water in 139cm column V;; air in 
water in 234cm column A; nitrogen in ethyl alcohol 139cm column oO: 
nitrogen in aqueous 0.5% PVA solution 139cm column @. 


- - -- Equation 6;- -—+—-+- Equation 7 combined with Equation 6. 
The following points and associated curves have been displaced by 
—4ems”'. 


Air bubbles rising in aqueous PVA solutions of various viscosities: 
54cP X; 115cP O; 216cP @ 


4,2. Discussion on rising velocities 


4.2. Discussion on rising velocities 


A force balance on a single bubble rising in a fluid may be written as: 
Net buoyancy force — Drag = Rate of change of bubble momentum. 


The last term must include the virtual mass of the bubble to account for 
that part of the momentum stored in the flow pattern of the surrounding 
liquid. Since a bubble of gas rising in water reaches steady state conditions 
very rapidly (i.e. within about one bubble diameter), [15], it is satisfactory 
to simplify the above expression by neglecting the unsteady state term. This 
simplification is also adequate for the cases of mass transfer considered in this 
paper since the rate of change in momentum due to mass transfer is small 
compared with the actual momentum associated with the bubble. Therefore 
the above equation may be written as 


4 
Fa = RTT (01 - 09) 9. (1) 


If the coefficient of drag (Cz) is defined on the basis of the projected area of 
an equivalent spherical bubble, then 


Fy 8 Gre Pg 4 
6x2 ae Gee ee (2) 
sa (are) <8 UP ( 


where Fr is a FROUDE number 
v2 
Fr= roe ees 
.{1 — 22). 
2 ( " ) IGTe 


Since £4 « 1, it can usually be neglected in the above expression. 

A dimensional analysis to find the dimensionless groups which may affect the 
velocity, providing the density and viscosity of the gas phase can be ignored, 
leads to a relationship between four such groups: 


f(Ca.%, We, =) = 0 (3) 


or 


f (Ca, We, Mo, = ) =0, (4) 


where Sf is the REYNOLDS number, Kt = rae We is the WEBER number, 


2 
We = arr and the parameter Mo represents the MORTON number, 


4 
an, 
Mo = 245. 


PLY 
As shown by HABERMANN and Morrow, [11], these groups are not sufficient 
to describe the behaviour of small bubbles completely. In the case of large 


4. RISING VELOCITY 


bubbles of spherical cap shape rising in liquids of low viscosity, the REYNOLDS, 
WEBER and MORTON numbers are not important. Form drag predominates 
and DAVIES and TAYLOR, [10], have shown experimentally that the velocity 
of these bubbles is given by 


U=-- g° Re, (5) 


where R, is the radius of curvature of the spherical cap. This velocity gives 
the correct pressure distribution in the immediate vicinity of the forward 
stagnation point, assuming that the flow in this region conforms to potential 
flow around a sphere. HABERMANN and MORTON obtained an almost 
identical result but expressed the velocity in terms of an equivalent radius 


U =1.02- Jg-re. (6) 


While it is shown that the results reported here are in agreement with 
previous investigations on spherical cap bubbles, it will be seen from Figure 1 
that the terminal velocity falls below that predicted by Equation 6. This 
discrepancy is attributable to wall effects since downward flow must occur 
between the bubble and the container walls. UNO and KINTNER, [20], carried 
out a number of experiments to investigate wall effect and developed the 
following expression, which agrees satisfactorily with the results of this 
investigation 
0.765 
Uy =109-U-(1- 5) (7) 
dt 

where U,, is the velocity in the presence of the wall and U is the velocity in 
an infinite medium. 

The retarding effect of viscosity above 54cP is not clearly understood. It 
is unlikely to be due to the proximity of the bubble to the walls since the 
effect is less the larger the bubble. As indicated in Figure 2, if the velocity 
is plotted versus the square root of the observed radius of curvature, all the 
points lie on the same curve. Thus Equation 5 appears to hold even with 
these viscous liquids. Since the equation was derived on the assumption of 
inviscid flow, it is unlikely to apply under more extreme conditions. The 
shape of the bubble obviously alters as the viscosity changes from 100 cP to 
200 cP, so that a relation which predicts velocities should take account of 
the viscosity of the fluid. DAVIDSON has also observed a decreasing radius of 
curvature and reduced velocity of spherical cap bubbles in viscous media, 
[7]. He attributes these observations to viscous forces within the body of the 
fluid, which result from the distortions of the stream lines around the body. 
There are other factors which merit consideration. The decreased radii of 
curvature in viscous solutions are likely to reduce the coefficient of form 
drag even though the overall coefficient of drag has obviously increased. 
The coefficient of skin friction drag may have increased due to the presence 
of surface active agents accumulating towards the rear of the bubble and 
producing an immobile interface. 


4.2. Discussion on rising velocities 


Rising velocity, cm/s 


1.0 1.4 1.8 2.2 2.6 3.0 3.4 3.8 4.2 4.6 
Radius of curvature, cm 


Figure 2: Velocities of bubbles as a function of radius of curvature. 
water A; 54cP V; 115cP x; 216cP @. broken line Equation 5 


The increased viscosity reduces the REYNOLDS number. For example, the 
REYNOLDS number based on the equivalent diameter’ for a 2cm diameter 
bubble rising in water is 6 x 10°. The corresponding values for the same size 
of bubble in 5.4% and 6.1% PVA solutions are 50 and 25 respectively. From 
the results of HABERMANN and Morrow, [11], the drag coefficient and the 
bubble shape appear markedly dependent on the REYNOLDS number when 
this is below 2 x 10° and moreover, the WEBER or the MORTON number 
is often an important additional parameter. This implies that in general 
both bubble shape and rising velocity are dependent on an interaction of 
hydrodynamics, surface tension and viscous forces. In the case of large 
bubbles in viscous solutions, surface tension forces are likely to be relatively 
small and hence it is not surprising that the viscosity influences the terminal 
velocity even though the bubble shape remains spherical cap. 


‘The REYNOLDS numbers based on the basal diameter are 11.5 x 10°, 95 and 36 for water for 
5.4% PVA solution and for 6.1% PVA solution respectively. 


5. BUBBLE SHAPE 


5. BUBBLE SHAPE 
5.1. Results 


The observed shapes and characteristics of bubbles rising in liquids are 
outlined in Table 2. The measured heights and basal radii for bubbles 
rising in water and PVA solutions are shown in Figure 3. The large bubbles 
were of the spherical cap shape observed by previous investigators, [10, 11]. 
Photographs showed that the curvature of the front surface of bubbles greater 
than 10cm” was very nearly spherical. On any particular bubble the distance 
from the front surface to the centre of curvature was within +4% of the 
mean radius. 

In low viscosity liquids the rear surface showed pronounced rippling but, 
with increasing PVA, these ripples were damped out. In solutions containing 
not less than 4.5% PVA, corresponding to viscosities greater than 54 cP, 
the rear surface became slightly concave. This could be due to a stagnation 
point produced by a trailing annular vortex. The maximum penetration at 
this point is 3mm, which corresponds to an upward velocity relative to the 
bubble of some 20cms ' to 24cems', which is similar to velocities in the 
trailing vortex which were observed visually. 

Large bubbles, greater than 15 cm”, rising in 6.1% PVA solution with a 
viscosity of 216cP, developed long, thin peripheral skirts as sketched in 
Figure 4. Measurements of the dimensions of the main body of the bubble 
indicated that only small amounts of gas were present in these skirts. The 
development of such skirts has not previously been reported. However, 
recent photographs of large bubbles rising in glycerine and oil, [1], suggest 
the presence of a trailing gas envelope. 


5.2. Discussion on bubble shape 


Dimensional analysis can also be used to relate bubble shape to other liquid 
properties. If the gas properties are not important and the only relevant 
liquid properties are density, viscosity and surface tension, the shape, which 
for a spherical cap bubble may conveniently be expressed in terms of a width 
to height ratio E, equal to 22, [15], will be related to the other properties by 


f (E,%, Mo, =) =0 (8) 


Since neither the surface tension nor the viscosity are likely to be important 
in the case of large bubbles, E should be constant for low values of a 
However, even though the results of the present tests show less scatter than 
those of DAviES and TAYLOR, [10], (Figure 3), there is still some doubt 
about the constancy of the value of £, which for an equivalent radius of less 
than 1cm appears to be influenced by liquid properties. 

With larger bubbles, equivalent radius 2cm, the ratio of width to height 
decreases with size as would be expected if the proximity of the walls were 


10 


5.2. Discussion on bubble shape 


+ Basal radius 
 einnaiiniammmmnemeemmmniee ane cal 


Linear dimensions height and basal radii, cm 


/3 


WO ‘liIpeJ jeseg pue 1yHlay SUOISUaLUIP JeOUIT 


0.5 1.0 1.5 2.0 2.5 
Radius of equivalent sphere, cm 


Figure 3: Heights and basal radii for rising bubbles as a function of equivalent 
radius. 
Lower graph (a) bubbles in water +; DAVIES and TAYLOR results [10] 
height m, basal radius @. 
Upper graph (b) bubbles in PVA solutions; 54cP O, 216cP x; lines as 
for water. 
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5. BUBBLE SHAPE 


(a) 
(b) 


Figure 4: Sketches showing relative shapes of 20cm* bubbles in (a) water and (b) 
216cP PVA solution 


influencing bubble shape. Figure 3 also indicates that as the concentration 
of PVA increases and the solutions become more viscous, the ratio of the 
width to the height again decreases. This effect, which implies a decreasing 
radius of curvature with increasing PVA, is most noticeable with smaller 
bubbles (r. * 1.7) in the 216 cP solution. In addition, as Table 2 shows, the 
size at which the spherical cap shape is adopted appears to decrease as the 
viscosity increases. 
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Liquid pupele os pe ie Remarks 
cm 
Tap water and distilled water 0.04—4.5 0.21-1.02 Ellipsoidal bubbles—marked oscillation. Helical rise. Transition to spherical 
cap at 4.5 cm°®. The transition from spherical cap to ellipsoidal is very 
noticeable during absorption of CO, from bubbles of suitable size range. 
4.5-40 1.02—2.1 Spherical cap bubbles. Rippling evident both on rear and front surfaces and 
on the relatively plane undersurface. 
Ethyl Alcohol <1 < 0.62. Ellipsoidal bubbles with marked oscillation. 
1-40 0.62—0.95 Spherical cap bubbles as for water. Rippling slightly more pronounced. 
0.5% PVA (4 = 2cP) <1 < 0.62 — Ellipsoidal bubbles with no oscillation. 
1-3.5 0.62—0.95 Ellipsoidal bubbles. Oscillations evident but less noticeable than with ethyl 
alcohol or water. 
3.5-40 0.95-2.1 Spherical cap bubbles-rectilinear rise. Rippling much less pronounced than 
with water. 
4.2% PVA (yu = 54cP) <1 < 0.62 — Ellipsoidal — no oscillation. 
1-3 0.62-0.89 ‘Transition region. Ellipsoidal with slight oscillation. 
3-40 0.89-2.1 Spherical cap bubble. Rectilinear motion. Apart from initial release rippling 
very slight. Lower surface slightly concave. 
5.4% PVA (u = 115 cP) <1 < 0.62 — Ellipsoidal. No surface instability. 
1-1.5 0.62-0.71 Slight surface instability. Gradual change of shape to spherical cap. 
1.5-40 0.71-2.1 Spherical cap. Rectilinear rise. No rippling. Lower surface concave. 
6.1% PVA (ps = 2.6 cP) 0.2-15 0.36-1.5 Spherical cap. Rectilinear rise. No rippling or oscillation. No shape transition 
in this range. 
15-40 1.5-2.1 Bubble tends to form thin skirts of trailing gas. Length of skirt increasing 


with bubble size, up to approx. 30cm long, 


Table 2: General description 


6. MASS TRANSFER 


6. MASS TRANSFER 
6.1. Results 


The experimentally determined rates of absorption have been interpreted in 
terms of a mass transfer coefficient, 6, defined as follows’, 


= -BA; + (Cy - Cp) . (9) 
A has been taken as the surface area of a spherical bubble with the same 


1 
volume as the real bubble and is equal to (36nV" Vee Since the gas contained 
only COxg, all the resistance to transfer was confined to the liquid. BAIRD 
and DAVIDSON, [2], have shown that, under isothermal conditions and with 
equilibrium at the interface, the rate of mass transfer may be given by the 
following equation, 


1 dV. Ps 
RO (Po: ty SB) = -asv 


3 
2 


n= -(Pco,:H-Cy) (10) 
1 
where S is the shape factor, in this case (367)? and H is the HENRY’s law 
constant. In the case of absorption from CO, bubbles into water free from 
COs, 
Poo, = Pp, Ch =0. (11) 

Furthermore, since the velocity is only slightly dependent on the bubble 
volume, it has been assumed constant at a mean value U. Hence 

dP = 

aro (12) 


Substitution in Equation 10 and integration with 6 constant leads to the 
following equation for the mass transfer coefficient in terms of other bubble 


parameters, 
ce hSH ZO 
VP8.vo-VP?.vt = a (Vt P2)4- VPP) A (13) 


This mass transfer coefficient may be related to the coefficient hf used by CALDERBANK and 
LOCHIEL, [3], by, 
* Ar 
B = hy ‘ “A? 


where Av is the total surface area of the bubble. The value of aa 


may be calculated from 
the ratio 22 From Figure 3(a), for bubbles between 1 cm and ne equivalent radius rising in 
water, this ratio may be shown to equal 3.90, which gives hj = 0.61- 8. For bubbles of the same 


size rising in 216 cP solution, the ratio 2 is 3.25 corresponding to a value of h7 = 0.68- 8. In 
spherical cap bubbles, with which this work is principally concerned, it is difficult to make a 
precise estimate of the proportion of mass transfer through the rear surface and consequently, it 
is not possible to separate completely the effect of changes in area from changes in mass transfer 
coefficient, which is the main reason for using h;. In addition, many of the bubbles exhibited 
a rippling surface which would make accurate estimation of the true area difficult. For this 
reason, the more common practice of using a coefficient based on the equivalent spherical area is 


preferred. 
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6.1. Results 


where the superscripts 0 and T refer to bubble properties at the time of 
release and at the top of the column respectively. 

In the case of desorption, Equation 9 is still applicable and Equation 10 
will also be valid if any transfer resistance in the gas phase is negligible 
compared with that in the liquid. However, Equation 11 no longer holds and 
the pressure of CO, in the bubble is given by 


RO PR V° 
Peo, = (n=) FE = [ry PEE) (14) 


For desorption, the liquid in the column was saturated by small bubbles of 
COzg rising from the bottom. This could have resulted in a gradient of CO, 
concentration along the column, of at most 20%, if no mixing had occurred 
during saturation or during each bubble experiment. Thus Cy could have 
varied with height according to the equation 


Ci, = Ps: H = Pe: HA, (15) 


and it is probably best to allow for this maximum variation. Substituting 
Equation 14 and Equation 15 into Equation 10 and integrating leads to 


3/ T wT\4 _ 3 0. o\4 _ PSH RO 
-(P3-v")-( VRE [) a0 


No allowance has been made in any of the calculations for the transfer of 
water vapour or nitrogen between the liquids and the bubbles. Approximate 
calculations show that these effects should be much less than the experimental 
errors arising from the measurement of velocities and volumes. It is believed 
that the derived mass transfer coefficients are accurate to about +20%. 
The mass transfer coefficients vs. the equivalent radius are shown in Figure 5 
and Figure 6. This radius is based on the arithmetic mean of the bubble 
volume at release and at the top of the column. Since the volume changes 
were relatively small and, as may be seen from the graphs, ( is not strongly 
dependent on r,, the use of the arithmetic average volume will be within the 
limits of experimental accuracy. It is possible that unsteady conditions at 
and immediately following bubble release may have influenced the results 
and have caused the mass transfer coefficients recorded in experiments with 
the 164cm column to be slightly higher than those from the 250 cm column 
(Figure 6). In view, however, of the comparatively good agreement with 
other investigators and the short duration of unsteady conditions, [15], such 
effects are not thought to be significant. 
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6. MASS TRANSFER 


Mass Transfer coefficient, cm/s 


0.046 


0.042 


0.038 


0.034 
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0.018 


oon 0.4 0.6 0.8 1.0 1.2 1.4 16 1.8 2.0 22 2.4 


Radius of equivalent sphere, cm 


Figure 5: Mass transfer coefficient for COg absorption in water and alcohol in 


164cm column: distilled water 19°C A; tap water 16°C @; ethyl alcohol 
O. Calculated curve for alcohol based on experimental results for water 
assuming h, VQ, broken line 
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6.1. Results 


Mass Transfer coefficient, cm/s 


0.014 


0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 
Radius of equivalent sphere, cm 


Figure 6: Mass transfer coefficients for CO, being absorbed by various liquids. 
Full lines: tap water absorbing; 164cm column, 16°C 0; 250cm column, 
15°C @. Broken lines: PVA solutions absorbing; 2cP x; 44cP A; 58cP 
A; 115cP V 
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7. DISCUSSION OF RESULTS 


7. DISCUSSION OF RESULTS 


7.1. Comparison with previous work 


The majority of the earlier investigations have been confined to smaller 
bubbles. However, experiments by LEONARD and HOUGHTON, [14], and 
BAIRD and DAvIDsOoN, [2], which have been shown in Figure 7, are in fair 
agreement with the present results where the ranges overlap. With the 
exception of results for a bubble of 1.25cm equivalent radius, CALDERBANK 
and LOCHIEL’s results” are higher than those obtained in the present in- 
vestigation. This could be at least partially attributable to their use of 
purer solution. Their experiments were confined to a 4in i.d. column; only 
specially purified distilled water was used and it was frequently changed. As 
shown later, surface active agents may considerably retard mass transfer 
even to large bubbles. 


7.2. Absorption and desorption 


In addition to the results in Figure 5 and Figure 6, the mean mass transfer 
coefficients obtained from measurements on CO, desorbing into nitrogen 
bubbles confirm the results from absorption and indicate that any resistance 
on the gas side of the interface is much less than on the liquid side. The 
desorption results were, however, less reproducible than the absorption 
experiments and this may be attributed to two causes. Firstly, some mass 
transfer occurs to the bubble while it is still in the cup so that the actual 
concentration difference is slightly less than the one used in the calculations. 
Secondly, a certain amount of mixing will be brought about by the bubble 
and some of the liquid in the bottom of the column will be carried up in the 
bubble wake. This will result in a concentration difference (Cy — C;) in this 
area, which is greater than the estimated one. Though these two errors act in 
opposite directions and may cancel out, they reasonably explain the scatter 
in the results. Another error may result from relating the CO» concentration 
in the liquid to the hydrostatic head (Equation 15) and neglecting the mixing 
effect of convection currents. 


’These have been recalculated to give a mass transfer coefficient based on the area of the equivalent 
sphere. 
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7.2. Absorption and desorption 
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Mass transfer coefficient, cm/s 
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Figure 7: Comparison of mass transfer coefficients (related to the area of the equiv- 
alent spherical bubble) for CO, being absorbed by water, calculated from 
the results of various investigations: 
this work shaded band; BAIRD and DAVIDSON [2] X; LEONARD and 
HOUGHTON [14] @; CALDERBANK and LOCHIEL [3] A, each point repre- 
senting a group of results which are close to one another; COPPOCK and 
MEIKLEJOHN [4| @, HAMMERTON and GARNER for the absorption of 
ethylene into water [12] O; Full lines A and B see text 
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7.3. The effect of surface active agents and viscosity 


Because PVA reduces the surface tension of water, it is not possible to 
separate completely the effects of viscosity and surface activity. Trace 
amounts of surface active materials present in tap water are presumably 
responsible for the reduction of the mass transfer coefficient below that for 
distilled water shown in Figure 5. Figure 6 indicates that small additions of 
PVA to distilled water reduce the mass transfer coefficient by approximately 
50%, although the viscosity is only raised to 2cP. Further additions of PVA 
cause a continuous increase in the mass transfer coefficient until 216 cP is 
reached. The curve for this viscosity (which is not shown in Figure 6) fell 
slightly below that for the 115 cP solution for bubbles up to 1.3 cm equivalent 
radius. Between 1.3cm and 1.8cm radius the mass transfer coefficients rose 
from 0.020cms ' to 0.027cms'’; they then fell to 0.026 cms for bubbles 
up to 2.1cm radius. 

The higher figures obtained in the 216 cP solution, for bubbles of 1.8cm and 
above are almost certainly attributable to the additional surface area caused 
by the formation of skirts (see Table 2). 

It is difficult to account for the effect of PVA concentration on mass transfer 
coefficient. One or more of the following factors may be important. 


1. Small quantities of PVA reduce rippling and this in turn reduces the 
mass transfer rate. A similar effect was observed by BAIRD and DAVID- 
SON, [2], using other surface active agents. However, this suppression of 
rippling fails to explain the increase of the mass transfer coefficients as 
the PVA concentration increases above 0.5 % for, with these solutions, 
rippling was damped even more effectively. 


2. The presence of an additional film resistance due to an accumulation 
of surface active agents at the interface. This latter effect has been 
noted previously for ripple-free flow around a solid object, [5]. How- 
ever, in the present series of experiments, the diffusion coefficients 
for COyg in the PVA solutions were measured using the spherical film 
absorption technique, [8]. In this method, solutions are passed over a 
spherical surface and the rate of absorption of CO, is measured. The 
diffusion coefficient is calculated on the basis of unsteady diffusion into 
an infinite liquid. The time for diffusion, which is that required for 
the fluid flow round the sphere, allows ample opportunity for surface 
active agents to concentrate at the gas-liquid interface. Their presence 
would be detected by a lowering of the apparent diffusion coefficient 
of COg. Since no such lowering was observed, as may be seen from 
Table 1, the presence of a film resistance to mass transfer is unlikely. 
However, measurements of the CO, diffusivity were not carried out at 
the concentrations of PVA (0.5% to 1%), which showed the maximum 
reduction in the mass transfer coefficient. 


3. The formation of an immobile surface film due to the accumulation of 
surface active agents. Such a film would cause an appreciable velocity 
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gradient in the liquid at the bubble surface and thus cause a reduced 
rate of transfer compared with that obtained when potential flow occurs 
around the front section. A surface film may only require the presence 
of a limited concentration of surface active agents, sufficient to form a 
monolayer. Normally the presence of larger amounts would probably 
have no additional effect. When, however, as in the case of PVA, 
increased quantities of surface active agents raised the viscosity, the 
shear forces on the bubble surface would also increase. These forces 
would tend to remove the immobile surface film and to restore relative 
movement between the gas and the liquid, at least over part of the 
frontal surface of the bubble. 


7.4. Theoretical models 


Theoretical models for mass transfer into spherical cap bubbles have been 
proposed by BAIRD and DAvIpsoN, [2], and by CALDERBANK and LOCHIEL, 
[3]. Both are based on unsteady state diffusion into the liquid flowing past 
the bubble. BAIRD and DAVIDSON relate their mass transfer coefficients to 
an area of equivalent volume sphere (as in this paper), whereas CALDERBANK 
and LOCHIEL relate theirs to the frontal area of the bubble. 

SZEKELY, [19], has measured absorption of the COz from the rear of a bubble 
held in a thin perspex dome. He reports that the rate of transfer per unit 
area through the rear surface is ; (+60 %) that through the front. As the 
ratio of the frontal surface area to the rear is approximately 2 (subtended 
half-angle 54°), it is evident that the transfer through the rear part of the 
bubble amounts to about 20% of the total. When this amount of transfer is 
deducted from the results obtained with distilled water, the mass transfer 
coefficients represented by line A in Figure 7 are obtained. These may be 
compared with the values predicted Equation 17 for frontal transfer by the 
equation of BAIRD and DAVIDSON, represented by line B in Figure 7. 


1 1 1 
B=0.975-de4-GD2-g!. (17) 


If the approach of LOCHTEL and CALDERBANK, [15], were adopted, Equa- 
tion 17 would be replaced by 


wil 1 i 
BH s de” +O Uay 


The use of this would reduce the values for the calculated mass transfer 
coefficients by approximately 3% and 9% for bubbles with r. equal to 1cm 
and 2cm respectively, rising in a 15cm dia. column. 

The measured values lie somewhat above the calculated values: the difference 
may well be caused by ripples in the bubble surfaces which would enhance 
the rates of mass transfer. 

The DAVIDSON and HARRISON model, [9], for mass transfer round spherical 
surfaces can be extended to account for mass transfer to a trailing skirt of 
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gas, for the case where this has the shape of a cylinder extending vertically 
downwards from the base of the bubble. The rate of mass transfer through 
the skirt, Aski¢4, and the concentration, C’, are given by 


® (dC 
Nskir =27Q- (ss) d® (18) 
‘ ®) dw W=0 
and 
C-C, 2 WIS —y 
=|]-—.- dn. 1 
CO d 


In the region of the skirt, V and ® are defined as follows: 


27 -W = 27R,-sin(a)-U,-y 
d®& 
— = R?-sin’(a)-U, 
Dy = R2-sin?(a)-U,-1 +O. 


1 is the length of the skirt and ®,; has the same value as given by DAVID- 
SON and HARRISON’s equation, C.10, [9]. However, the definition of 2 is 
slightly different from that used by DAVIDSON and HARRISON, [9], because 
in the region of the skirt the stream lines are assumed to be parallel to 
the vertical axis. The slight change in the flow direction resulting from the 
approximation of the skirt to a vertical cylinder has been assumed to have 
a negligible effect on the concentration gradient. The surface velocity U, 
has been taken as constant over the length of the skirt and equal to J2gb 
or [29R. -y1l- cos(a) | where 0 represents the height of the bubble.“ With 


these modifications the overall mass transfer rate becomes 


= —4-Vn-D-8-(C)-G) . (20) 


The quantity of gas in the skirt is small and may be neglected when relating 
the radius of curvature of the spherical cap and the subtended half-angle 
to the bubble volume or equivalent diameter. On this basis, Equation 20 
can still be evaluated for the case of bubbles with skirts, to give the mass 
transfer coefficient based on the area of the equivalent sphere as 


=H 1 
B=K-det-V9- G4 (21) 
where 
Jims: (+) + 1.507 — 0.2667 » 2 — 0.6667+ ms? |‘ 
K=1,78.)| ——" (22) 
(m? + q)4 


“For bubbles subtending a 50° half-angle, the surface Velocity calculated from 
[29R. -V1l- cos(a) | is 0.84 + ./gR.. However, in the case of long skirts, this value will be 
an overestimate as the surface velocity must approach the terminal velocity given by Equation 5, 
ignoring wall effects. Since the mass transfer rate depends on the square root of the surface 
velocity, the estimate for the mass transfer coefficient over the skirt will not be more than 10% 
too high. 


22 


7.5. The effect of solute diffusivity 


and 


m = 1-cos(a) 
s=1+cos(a) 


q = 2+ cos(a) 


and a is the half angle subtended by the bubble at the centre of curvature. 
For the 50° subtended half-angle found with water, the coefficient is given by 


K = \/40-( ae) +095 (23) 


For a 65° subtended half-angle, which is applicable to bubbles in the more 
viscous solutions, the coefficient becomes 


k= yfo-(£) +00 2 


Measurements of skirt lengths in the 216 cP solution gave values for + of 
about 1.4 for 20mL bubbles and 1.8 for 30mL bubbles. For a subtended 
half-angle of 65°, the mass transfer coefficients for such bubbles should be 2 
to 2.5 times those for bubbles without skirts. This calculation neglects mass 
transfer from the inside surface of the skirt which would further increase this 
figure. 

Experimental measurements showed that mass transfer coefficients for bub- 
bles with skirts were only about 1.3 times those without. These provide 
further evidence that mass transfer at the interface is retarded by a film which 
either offers physical resistance to mass transfer or changes the hydrodynamic 
pattern so as to produce low velocities at the gas-liquid interface. 


7.5. The effect of solute diffusivity 


The models of HIGBIE, DANCKWERTS, BAIRD and DAVIDSON and CALDER- 
BANK and LOCHIEL all indicate the mass transfer coefficient to be propor- 
tional to the square root of the diffusivity. The results of HAMMERTON and 
GARNER for small bubbles of hydrogen (d, = 0.7cm) dissolving in water 
confirm this relation. 

If the results with distilled water (Prsec =16x10° cm s) are taken 


as a basis, then the assumption of 8 « VJ may be used to predict the 
values of mass transfer coefficients for bubbles of COg rising in ethyl alcohol 
(Airc =2.7x 10° cm” ae! As may be seen from Figure 5, the experi- 
mental results are similar to the calculated results for bubbles of 0.8cm to 
1.0cm equivalent radius. At larger radii the experimental results increase 
above the predicted values. 
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A. ALTERNATIVE EVALUATION OF Figure 2 


A. ALTERNATIVE EVALUATION OF 
Figure 2 


If Figure 2 is redrawn (see Figure 8) and subjected to regression analysis, 
the following relationship will be obtained with a coefficient of determination 
Of 0.963: 


Ces. 


This may be approximated by the relation 


2 
U=22.4-R2. 


Bubble rising velocity, cm s* 


1,0 1,5 2,0 2,5 3,0 3,5 4,0 4,5 
Radius of curvature, cm 


Figure 8: Rising velocity of bubbles in water plotted against radius of curvature 
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LIST OF SYMBOLS 


A surface area of a sphere of the same volume as the 
bubble 
basal radius of the bubble 
b height of bubble 
C concentration of CO» 
Ca coefficient of drag 
J diffusion coefficient of CO2 in appropriate solution 
d. diameter of sphere of equivalent volume to bubble 
d, diameter of column 
E ~ width to height ratio of spherical cap bubble, 2a 
Fy drag force 
Fr FROUDE number, a 
g acceleration due to gravity 
H ~~ HEnrRy’s Law constant 
K coefficient 
l length of skirt from base of bubble 
Mo MORTON number, a 
n moles of CO» in bubble 
n moles of CO, transferred to the bubble per unit time 
P pressure 
& gas constant 
R, radius of curvature 
9 REYNOLDS number, meee 


re equivalent radius (i.e. radius of sphere of equivalent 
volume to bubble) 


wle 


shape factor, (367) 
time 


S 
t 
U terminal velocity of bubble in an infinite medium 
U mean rising velocity of bubble 

Us 


velocity of liquid at the gas-liquid interface relative 
to the bubble 

U,, velocity of a bubble influenced by proximity of 
column walls 


V volume of bubble 


Urd 
We WEBER number, = 
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LIST OF SYMBOLS 


y distance along outward normal from bubble or skirt 
surface 

Zz vertical co-ordinate 

Qa half-angle subtended by bubble at centre of curvature 

B mass transfer coefficient 

y surface tension 

m dynamic viscosity of liquid 

p density 

0 temperature 

b bulk of solution 

B bubble 

gas 

i interface 


hs hydrostatic 
1 liquid 


at release 


at top of column 
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